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We present the ingredients for determining tt'^-k" matrix elements via the combination of lattice QCD 

and chiral perturbation theory (^PT). By simulating these matrix elements at unphysical kinematics, it is possible 
to determine all the low-energy constants (LECs) for constructing the physical K'^ tt^tt" amplitudes at next- 
to-leading order (NLO) in the chiral expansion. In this work, the one-loop chiral corrections are calculated for 
arbitrary meson four-momenta, in both xPT and quenched xPT (qxPT), and the finite- volume effects are studied. 



1. INTRODUCTION 

The need for a high-precision prediction for 
K TTTT amplitudes is underlined by the re- 
cent experimental measurement of Re(e'/e) and 
the long-standing puzzle, the AI =1/2 rule. Al- 
though the finite-volume (FV) techniques devel- 
oped in Refs. |^,^|J§1 can ultimately enable an 
accurate calculation oi K tttt decay rates, the 
most practical approach to the numerical calcu- 
lation of these decay rates remains the combina- 
tion of lattice QCD and (quenched and partially 
quenched) xPT. Apart from a calculation for the 
CP-conserving, AI = 3/2, K ^ tttt decay in 
Ref. ||l|, all the numerical studies hitherto fol- 
low a strategy proposed in Ref. which only 
allows the determination of these amplitudes at 
leading-order (LO) in the chiral expansion]^ Be- 
cause of the large kaon mass and the presence of 
final state interactions, non-LO corrections in this 
expansion are significant. In a recent work 
we have proposed to perform lattice simulations 
at unphysical kinematics over a range of meson 
masses and momenta, from which we can deter- 
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^In Ref. pi], the decay amplitude is also obtained at the 
precision of leading-order in the chiral expansion. 



mine all the necessary LECs for constructing the 
physical matrix element {'k'^tt^\0'^^^^\K^) at 
NLO in the chiral expansion. We have suggested 
a specific unphysical kinematics^ the SPQR kine- 
matics as explained in detail in Refs. which 
enables such a procedure, and have studied the 
following AS = 1 operators (a, /3 are colour in- 
dices) 

Qi = iSada)L{Uf3Ui3 - d/3d/3)L 
+ {SaUa)L{updf3)L, 

Qt = ■^{Sada)L ^ eg{qpq(3)R, (1) 
q—u,d.s 

2 ^ 

Qs = -j^{Sadp)L ^ eq{qf3qa)R, 
q—u,d,s 

where is the electric charge of q and {ipiip2)L,R 
means ^^i7^(l T 75)V'2- 

2. FINITE- VOLUME EFFECTS 

In Ref. 1^, we investigate the FV corrections, 
power-like in 1/L, which arise from replacing 
sums by integrals in the one-loop calculation that 
involves the diagrams in Fig. 1 ^. It can be shown 

''Another choice is considered in Ref. flq]. 

^ Such a calculation for (7r+7r'^|Q4|_ft'+y at two particular 

kinematics, Mj^ = M-,^ and Mx = 2Mt^ with all mesons 
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that a diagram which does not have an imaginary 
part in Minkowski space will only have FV correc- 
tions exponential in L, therefore only diagram (c) 
contributes to the 1/L" corrections. Because the 
two-pion final state has 1 = 2, this diagram only 
contains four-quark intermediate states and there 
are no disconnected quark-loops in the quark-flow 
picture. For the same reason, it does not receive 
contributions from the 77' propagator. Hence the 
effects are identical in xPT and quenched 
XPT (qxPT) for K+ ^ 



tt" at this order, and 



only at this order. This is not true for A/ =1/2 
decay amplitudes PpOl. 



(a) 



(b) 



(c) 




Figure 1. One-loop diagrams for K'^ — > tt+tt'' 
amplitudes. The grey circles (squares) are weak 
(strong) vertices. The diagrams for wavefunction 
renormalisation are not shown here. 



in the center-of-mass 
tt+tt" amplitudes, these 



We have found that 
frame, for all the 
one-loop 1 /L" corrections are independent of the 
weak operators and can be removed by a univer- 
sal factor derived first by Lellouch and Liischer 
in Ref. This modifies the conclusion of Refs. 
]ll| , |T^ , in which a FV effect resulting in the shift 
of the two-pion total energy in the argument of 
the tree-level amplitudes is interpreted as a gen- 
uine 1/L" correction to the matrix elements, and 



therefore the FV effects in 
found to depend on Mk- 



{TT+Tr°\Q4K+) are 



at rest, have been performed in Refs. 



We are currently investigating the FV effects 
of these amplitudes in a moving frame. The 
Lellouch-Liischer factor has not yet been derived 
for this, while the modification of Liischer's quan- 
tisation condition ||l^,|lj| relating the infinite- 
volume TTTT scattering phase to the FV two-pion 
energy spectrum, due to the moving frame was 
obtained in Ref. ||l^. As a by-product of our 
work, we verify that the energy shift obtained in 
one-loop perturbation theory in a moving frame 
agrees with the expansion of the quantization 
condition in Ref. jl^ to the same order. 

3. ONE-LOOP CHIRAL CORRECTIONS 
IN INFINITE VOLUME 

We evaluate the one-loop correction by us- 
ing dimensional regularisation and subtracting 
log(47r) - 7£; 1 -h 2/(4 - d). The lowest-order 
amplitudes are all proportional to 1//^, where / 
is the light pseudoscalar meson decay constant 
in the chiral limit. At NLO, we choose to ex- 
press 1//^ in terms of l/(/^/;f). This factor 
fully absorbs the dependence upon the Gasser- 
Leutwyler LECs L4 and L5 introduced via wave- 
function renormalisation. 

In Ref. the one-loop diagrams have been 
calculated for arbitrary external meson four- 
momenta in both xPT and qxPT. The results 
are lengthy and are presented on a web site [l^ . 
In Fig. 2, we show an example of these results 
for {tt~^tt'^\07^s\K^). These plots are the ratios 
between the one-loop corrections and the lowest- 
order matrix elements with both final-state pions 
at rest. Fig. 2a is the result in xPT and Fig. 
2b is that in qxPT. This figure suggests that in 
a quenched numerical calculation of these matrix 
elements, it is not appropriate in general to per- 
form chiral extrapolation using unquenched xPT 
results JT?! . This is confirmed by numerical data 



4. CONCLUSIONS 

We have made theoretical progress towards the 
calculation of K tttt decay amplitudes via the 
combination of lattice QCD and xPT. We find 
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Figure 2. Ratio between the one-loop correc- 
tion, at the renormalisation scale 0.7 GeV, and 
the lowest-order result for {tt^tt'^IOy.sIK'^) in (a) 
XPT and (b) qxPT. The x axis is and the y 
axis is Mk- In (b), the coupling accompanying 
the kinetic term of the rj' propagator is set to zero, 
and the rj' mass is taken to be Mq — 0.5 GeV. The 
one-loop results are not very sensitive to these pa- 
rameters. The singular behaviour along the line 
= \/2Mk in (b) is due to the fact that when 
performing the qxPT calculation, we use a ba- 
sis in which the pseudo-Goldstone states are qq' 
mesons, where q and q' are u, d and s, and the ss 
meson becomes tachyonic when M-^ > V^Mk. 



it feasible to determine all the LECs for con- 
structing /=2(7r7r|e'^'5=i|ii:) at NLO in the chi- 
ral expansion. A quenched numerical study is in 
progress As for the A/ = 1/2 channel, we 

find the situation to be considerably more com- 
phcated 
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